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The  calibration  data  compiled  and  analyzed  in  this 
report  were  taken  from  Advanced  Development  (AD)  Program 
Inertial  Measurement  Units  (IMU)  during  the  March  1976  to 
October  1978  time  frame.  Since  that  time  some  requirements 
and  specifications  given  in  the  body  of  the  report  and 
various  tables  and  charts  have  been  changed  for  the 
Engineering  Development  (ED)  Program  based  on  findings 
during  AD  and  pre-ED  work.  The  criteria  and  discussions 
used  in  this  report  are  those  relative  to  AD. 

During  the  pre-ED  time  period  three  IMUs  (serial 
numbers  007,  Oil,  and  012)  were  re-worked  to  incorporate 
many  of  the  improvements  identified  during  AD.  As 
anticipated  the  resulting  calibration  parameter  behavior  is 
much  improved.  A  future  report  will  document  the  changes 
and  analyze  the  performance  of  the  modified  units. 
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I.  INTRODUCTION  AND  OBJECT 

The  Pershing  PII  Inertial  Measurement  Unit  (IMU)  is  a 
four  gimbal,  all  attitude  system  incorporating  two  two- 
degree-of-f reedom  dry  tuned  gyros,  one  dual  axis  and  one 
single  axis  force  rebalance  accelerometer. 

IMU  performance  is  largely  dependent  on  gyro  and 
accelerometer  parameters  which  must  either  exhibit  a  high 
degree  of  absolute  accuracy  over  the  usable  life  of  the 
IMU,  or  which  must  be  amenable  to  a  periodic  calibration. 

The  PII  IMU  is  a  calibrated  system  and  as  such  abso¬ 
lute  accuracy  of  gyro  and  accelerometer  parameters  is  unim¬ 
portant.  Instrument  channel  errors  are  measured  in  the 
calibration  procedure  and  are  software  compensated  so  that 
a  high  degree  of  parameter  stability  is  required  rather 
than  absolute  accuracy. 

A  significant  advantage  of  a  calibrated  system  is 
cost  savings  resulting  from  relaxation  of  the  requirements 
to  design,  fabricate  and  test  for  high  absolute  accuracy  at 
both  the  component  and  system  level.  Instead,  nominal 
parameter  values  are  built  in  and  are  subsequently  deter¬ 
mined  very  accurately  via  the  calibration  procedure,  allow¬ 
ing  software  compensation  to  be  applied  during  ground 
alignment  and  flight  modes.  Frequency  of  calibration  is 
determined  by  the  stability  of  the  measured  parameters. 

Some  of  the  advantage  of  a  calibrated  system  is  therefore 
lost  if  stability  is  such  that  frequent  and  lengthy  cali¬ 
brations  are  required. 

The  object  of  this  study  is  to  compile  calibration 
data  from  all  Engineering  Model  (EM)  and  Prototype  Model 
(PM)  IMUs  for: 

(1)  Identification  of  least/most  stable  parameters 

(2)  Application  of  acceptance  criteria  to  all  runs- 

(3)  Identification  of  data  trends  which  may  be  useful 
to  predict  the  time  interval  between  calibrations. 

The  above  information  should  be  helpful  in  determin¬ 
ing  if  the  currently  proposed  calibration  frequency  of  once 
per  six  months  is  realistic  or  if  the  time  interval  should 
be  increased  or  decreased.  The  data  should  provide  visi¬ 
bility  of  parameters  which  dictate  a  minimum  time  between 
calibrations.  Application  of  acceptance  criteria  to  all 


calibration  results  should  yield  further  insight  into 
parameter  behavior  and  overall  system  performance.  An 
analysis  of  data  trends  should  also  aid  in  corroborating 
parameter  stability  and  time  between  calibrations. 

II.  DISCUSSION  OF  CALIBRATION  PROCEDURE  AND  OUTPUT  DATA 

Calibration  of  the  PI I  IMU  is  an  automatic  proce¬ 
dure^-  performed  under  computer  control.  The  instrument 
cluster  (azimuth  gimbal)  is  torqued  to  various  positions  to 
orient  accelerometers  relative  to  earth's  g  field  and  to 
orient  gyros  relative  to  earth's  rotational  field.  Data, 
in  the  form  of  accelerometer  and  gyro  torquing  pulses,  are 
acquired  in  a  total  of  13  positions.  These  data  are 
utilized  to  compute  34  gyro  and  accelerometer  parameters 
which  are  stored  in  the  IMU  Serial  Core  Memory  (SCM)  for 
later  use  in  the  ground  alignment  and  flight  modes. 
Currently,  approximately  4  1/2  hours  are  required  to 
perform  a  calibration. 

The  34  calibrated  parameters  are  shown  in  Table  1 
along  with  definitions  and  units.  Parameters  with  the 
letter  H  as  a  suffix  are  used  in  the  ground  alignment 
(gyrocompass)  mode  only.  DELTA  YX,  DELTA  ZY  and  DELTA  ZX 
are  factory/ depot  calibrated  parameters  and  are  measured  in 
subsequent  calibrations  for  goodness  checks  only.  These 
three  parameters  along  with  two  other  factory/depot  deter¬ 
mined  parameters,  HSX  and  HSY,  are  stored  in  a  protected 
portion  of  the  SCM  and  theoretically  are  never  updated 
except  at  the  factory/depot.  HSX  and  HSY  are  required  for 
accurate  gyrocompass ing  and  are  determined  by  a  special 
purpose  factory /depot  test  apart  from  the  standard  calibra¬ 
tion  procedure.  Factory/depot  determination  of  more  accu¬ 
rate  values  of  KIXH  and  KIYH  for  gyrocompass  purposes  has 
been  proposed.  This  also  requires  a  special  purpose  fac¬ 
tory/depot  test.  The  resulting  values  would  be  stored  in 
the  protected  portion  of  the  SCM  and  would  be  updated  only 
at  the  factory/depot. 

A  measure  of  parameter  stability  is  obtained  by 
taking  the  difference  between  the  value  obtained  from  a  new 
calibration  and  the  value  obtained  from  the  original 
reference  calibration.  The  34  parameter  differences  are 
defined  by  adding  the  letter  E  as  a  suffix  to  each  of  the 
measured  parameters,  e.g.,  DFZE  *  DFZ  (NEW)  -  DFZ  (REF). 
Table  2  shows  the  proposed  one  sigma  limits  placed  on  each 

Yl  The  Singer  Company,  Kearfott  Division*  Document 

Y256A337,  "Alignment  Gyrocompass  and  IMU  Parameter 
Compensation. " 


parameter  difference.  The  specification  requires  that  the 
IMU  meet  performance  criteria  with  a  calibration  frequency 
of  once  every  six  months.  It  is  therefore  of  prime  impor¬ 
tance  that  parameters  remain  stable  over  this  interval  of 
time  as  a  minimum. 

The  acceptance  test  criteria  utilized  during  the 
Advanced  Levelopment  (AD)  program  allowed  some  degradation 
of  the  limits  indicated  in  Table  2.  For  acceptance,  all 
parameter  values  were  required  to  lie  within  +3  sigma 
limits.  If  one  or  more  of  the  flight  parameters  exceeded 
+1  sigma  but  did  not  lie  outside  +3  sigma,  the  Root  Mean 
Square  (RMS)  value  of  all  flight  parameters,  each  normal¬ 
ized  to  its  one  sigma  value,  could  not  exceed  a  value  of 
1.3  for  acceptance.  Some  flight  parameters  could  therefore 
exceed  the  +1  sigma  limits  as  could  all  of  the  ground 
alignment  parameters  but  in  no  case  was  a  value  greater 
than  +3  sigma  acceptable. 

III.  DISCUSSION  OF  RESULTS 

Calibration  data  considered  in  this  study  were  taken 
from  a  total  of  13  IMUs  (4  EMs  and  9  PMs)  over  a  period  of 
30  months.  Calibration  sites  were  the  Guidance  and  Control 
Directorate,  Redstone  Arsenal,  Alabama  and  Martin  Marietta 
Corporation,  Orlando,  Florida.  The  reference  calibration 
was  ordinarily  performed  at  the  Singer  Company,  Kearfott 
Division,  Little  Falls,  New  Jersey,  however,  some  units 
required  re-referencing  at  the  other  sites,  usually  due  to 
electronic  card  replacement. 

Table  3  gives  a  tabulation  of  number  of  calibrations 
and  time  frames  involved  for  each  IMU. 

IMUs  with  serial  numbers  001  through  004  are  EMs  and 
the  remainder  are  PMs.  Listing  of  an  IMU  more  than  once 
indicates  that  more  than  one  reference  calibration  was 
utilized. 

Several  of  the  IMUs  as  indicated  in  Table  3,  had  very 
few  calibrations  performed,  particularly  those  that  were 
expended  in  AD  flights.  Only  those  units  with  six  or  more 
calibrations  over  a  period  of  six  or  more  months  were 
subjected  to  the  trend  analysis  presented  in  the  appendix. 
Otherwise,  a  total  of  117  calibration  runs,  which  comprises 
the  majority  of  runs  made  at  Redstone  and  Martin,  were 
considered  for  this  study. 
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Results  £rom  the  10  calibration  runs  made  with  IMU 
serial  number  005  during  11.4  months  period  are  presented 
in  Table  4,  as  a  data  sample.  The  data  have  been  normal- 
ized  by  dividing  the  measured  values  of  a  given  parameter 
by  the  allowable  one  sigma  value  given  in  Table  2. 

Table  5  shows  the  number  of  times  parameter  differ¬ 
ences  fell  In  the  indicated  sigma  ranges  for  each  IMU.  For 
example,  the  DFZE  parameter  from  IMU  serial  number  001  was 
within  the  prescribed  +1  sigma  limits  18  out  of  the  26 
calibrations.  The  parameter  fell  between  +1  and  +2  sigma 
four  times  and  between  +2  and  +3  sigma  four  times.  None  of 
IMU  serial  001  calibration  runs  yielded  a  DFZE  greater  than 
+3  sigma. 

Percentage  of  calibration  runs  falling  in  the  various 
sigma  ranges  for  each  parameter  for  all  IMUs  is  shown  in 
the  last  row  of  Table  5.  A  study  of  this  table  on  an  IMU 
basis  over  all  parameters  or  on  a  parameter  basis  over  all 
IMUs  givas  the  indication  that  certain  gyro  parameters 
exhibit  the  least  stability  whereas  certain  accelerometer 
parameters  are  the  most  stable.  These  parameters  are  dis¬ 
played  in  Figures  1,  2,  and  3  which  are  derived  from  Table 
5. 


Figure  1  depicts  the  percentage  of  parameters  from 
all  calibrations  which  fell  into  the  0+1  sigma  range. 
Parameters  in  the  left-most  columns  of  Figure  I  are  the 
least  stable  with  DIXE  and  DIYE  occurring  in  the  0+1  sigma 
range  for  only  43  percent  of  the  total  of  117  measurements 
for  each. 

The  range  is  expanded  to  0+2  sigma  in  Figure  2  which 
displays  the  same  12  parameters  In  the  two  left-most  col¬ 
umns  as  are  indicated  in  the  three  left-most  columns  of 
Figure  1. 

Figure  3,  in  which  the  range  is  increased  to  0+3 
sigma,  indicates  that  six  (left-most  column)  of  the  T2  ■ 
parameters  identified  in  Figures  1  and  2  are  the  very  worst 
performers.  The  six  are  all  gyro  parameters,  the  defini¬ 
tions  and  units  of  which  are  given  in  Table  2. 

The  most  stable  parameters  are  the  nine  appearing  in 
the  right-most  column  of  Figure  1.  They  are  all  acceler¬ 
ometer  parameters.  Definitions  and  units  for  these  are 
given  in  Table  2  also. 
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Placement  of  the  parameters  in  the  vertical  order 
shown  in  Figures  1  through  3  has  no  significance. 

Data  in  Table  5  are  further  condensed  in  Table  6 
which  provides  a  summary  of  frequency  of  occurrence  of 
difference  parameters  within  the  indicated  sigma  ranges  for 
all  IMUs  and  all  calibrations.  All  parameters  from  all 
IMUs  fall  within  the  prescribed  one  sigma  limits  for  only 
75  percent  of  the  total  number  of  data  points,  however,  91 
percent  lie  within  0+2  sigma  bounds  and  96  percent  lie 
within  0+3  sigma  bounds. 

Table  7  presents  results  of  the  application  of  AD 
acceptance  test  criteria  to  all  117  calibration  runs 
considered  in  this  study.  Fifty-five  of  117  runs  or  47.0 
percent  passed  acceptance  criteria,  whereas  11  runs  or  9.4 
percent  failed  all  criteria.  One  or  more  flight  parameters 
greater  than  +3  sigma  in  combination  with  larger  than 
acceptable  flTght  parameter  RMS  values,  accounted  for  the 
highest  combination  failure  category  of  24  runs  or  20.5 
percent. 

Table  8  lists  RMS  values  of  the  20  flight  parameters 
for  each  calibration  of  each  IMU.  Sixty-five  of  117  runs 
or  55.6  percent  passed  the  RMS  criterion  of  1.3  sigma.  The 
remaining  runs,  with  one  exception,  yielded  one  or  more 
parameters  greater  than  +3  sigma  and  thus  failed.  The 
single  exception  is  run  Number  9  from  serial  number  001 
which  failed  to  meet  the  RMS  criterion  of  1.3  sigma  even 
though  no  parameter  greater  than  +3  sigma  was  measured. 
Thirteen  runs  or  11.1  percent  met  the  RMS  criteria  even 
though  one  or  more  parameters  were  greater  than  +3  sigma. 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  taken  the  initial  overview  of  the 
majority  of  calibration  data  from  all  AD  IMUs. 

Identification  of  parameters  which  most  often  exceed¬ 
ed  three  sigma  limits  has  been  accomplished.  These  are 
DFX,  DIX,  diy,  DIZ,  KTX  and  KTZ  which  are  all  gyro  para¬ 
meters  .  Most  stable  parameters  are  KIZH,  KOY,  KOZ,  KOXH, 
KOYH,  KOZH,  DELTA  YX,  DELTA  ZX  and  DELTA  ZY  which  are 
associated  with  the  accelerometers. 

Application  of  AD  acceptance  criteria  to  all  runs 
shows  that  only  47.0  percent  passed  the  three  sigma 
criterion  and  55.6  percent  passed  the  1.3  sigma  flight 
criterion. 
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The  Appendix  addresses  the  question  of  trends  i>y 
using  the  linear  least  squares  regression  method.  The  three 
sigma  criterion  was  used  to  establish  a  predicted  worst 
case  time  between  calibrations  and  in  very  few  cases  is 
this  time  less  than  the  AD  program  required  time  of  six 
months.  However,  many  of  the  parameters  would  fail  the  six 
months  time  criterion  when  one  sigma  limits  are  imposed. 

The  plotted  data  in  *;he  Appendix  and  the  correspond¬ 
ing  tables  give  some  visibility  on  trending  parameters. 
Parameters  with  correlation  coefficients  near  unity  are 
trending  suspects.  Gyro  parameters  KTX,  KTY  and  KTZ  and 
accelerometer  parameters  Kix,  KIXH,  KIY  and  KIYH  show  the 
greatest  propensity  toward  trending.  Conversely,  those 
parameters  with  correlation  coefficients  near  zero  are  more 
random  in  nature. 

The  general  conclusion  is  that  many  of  the  parameters 
need  to  be  better  behaved.  There  appears  to  be  some 
problem  associated  with  test  locations  as  indicated  by  the 
plotted  data  in  the  Appendix  from  ZMU  serial  number  005, 
i.e.,  note  data  points  taken  at  Martin  (M)  compared  to 
those  taken  at  the  Guidance  and  Control  Directorate  (A). 

It  is  recommended  that  close  surveillance  of  the 
calibration  parameters  be  continued  into  the  Engineering 
Development  (ED)  program  as  component  and  system  improve¬ 
ments  are  incorporated.  Those  parameters  which  are  not 
well  behaved  should  be  reviewed  for  their  effects  on 
alignment  and  flight  accuracy,  and  methods  for  improving 
performance  should  be  defined  and  implemented. 

A  factor  which  requires  further  investigation  is 
minimization  of  time  required  for  performing  calibration. 
The  current  routine  has  not  been  time-optimized  for  field 
usage.  Elimination  of  certain  parameters  required  only  in 
a  factory  calibration  and  optimizing  the  sequence  of  posi¬ 
tions  to  avoid  backtracking  between  positions  are  two  pos¬ 
sibilities  for  decreasing  time. 

Additionally,  evaluation  of  the  calibration  operation 
under  field-type  conditions  will  be  necessary  to  determine 
the  impact  of  environmental  factors  and  unknown  initial 
heading  on  calibration  time. 
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PARAMETER 
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Figure  1.  Percent  of  parameters  in  the  0+1  sigma  range. 
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TABLE  1.  CALIBRATED  PARAMETERS 


Parameter 

KTX 

KTY 

KTZ 

DFZ 

DSZ 

DIZ 

KIY 

KIYH 

KIX 

KIXH 

DOZ 

DFX 

KIZ 

KIZH 

DFY 

DSX 

DSY 

DIX 

DIY 

DELTA  YX 

DELTA  ZY 

KOZ 

KSZ 

KOZH 

KSZH 

KOX 

KSX 

KOXH 

KSXH 


X  Gyro  Torquer  Scale  Factor 
(deg/hr  per  deg/hr) 

Y  Gyro  Torquer  Scale  Factor 
(deg/hr  per  deg/hr) 

Z  Gyro  Torquer  Scale  Factor 
(deg/hr  per  deg/hr) 

Fixed  Drift  of  Z  Gyro  (deg/hr) 

G  Sensitive  Drift  of  Z  Gyro  (Spin  Axis) 
(deg/hr  per  ft/sec2) 

G  Sensitive  Drift  of  Z  Gyro  (Input  Axis) 
(deg/hr  per  ft/sec2) 

Y  Accelerometer  Scale  Factor  (Flight)  (ug/g) 

Y  Accelerometer  Scale  Factor  (Ground)  (ug/g) 
X  Accelerometer  Scale  Factor  (Flight)  (ug/g) 
X  Accelerometer  Scale  Factor  (Ground)  (ug/g) 
G  Sensitive  Drift  of  2  Gyro  (Output  Axis) 
(deg/hr  per  ft/sec2) 

Fixed  Drift  of  X  Gyro  (deg/hr) 

Z  Accelerometer  Scale  Factor  (Flight) (ug/g) 

Z  Accelerometer  Scale  Factor  (Ground) (ug/g) 
Fixed  Drift  of  Y  Gyro  (deg/hr) 

G  Sensitive  Drift  of  X  Gyro  (Spin  Axis) 
(deg/hr  per  ft/sec2) 

G  Sensitive  Drift  of  Y  Gyro  (Spin  Axis) 
(deg/hr  per  ft/sec2) 

G  Sensitive  Drift  of  X  Gyro  (Input  Axis) 
(deg/hr  per  ft/sec2) 

G  Sensitive  Drift  of  Y  Gyro  (Input  Axis) 
(deg/hr  per  ft/sec2) 

Non-Orthogonality  Between  Y  &  X  Acceler¬ 
ometers  (rad) 

Non-Orthogonality  Between  Z  &  X  Acceler¬ 
ometers  (rad)  2 

Z  Accelerometer  Bias  (Flight)  (ft/sec  ) 

Z  Accelerometer  Scale  Factor  Asymmetry 
(Flight)  (ug/g)  2 

Z  Accelerometer  Bias  (Ground)  (ft/sec  ) 

Z  Accelerometer  Scale  Factor  Asymmetry 
(Ground)  (ug/g)  2 

X  Accelerometer  Bias  (Flight)  (ft/sec  ) 

X  Accelerometer  Scale  Factor  Asymmetry 
(Flight)  (ug/g)  2 

X  Accelerometer  Bias  (Ground)  (ft/sec  ) 

X  Accelerometer  Scale  Factor  Asymmetry 
(Ground)  (ug/g) 


TABLE  1.  (Concluded) 


Non-Orthogonality  Between  Z  &  Y  Acceler¬ 
ometers  (rad)  2 

Y  Accelerometer  Bias  (Plight)  (ft/sec  ) 

Y  Accelerometer  Scale  Factor  Asymmetry 

(Flight)  (ug/g)  2 

Y  Accelerometer  Bias  (Ground)  (ft/sec  ) 

Y  Accelerometer  Scale  Factor  Asymmetry 
(Ground)  (ug/g) 


TABLE  2.  DIFFERENCE  PARAMETER  LIMITS 


Difference 

Parameter  +lo  Max  Criteria  +3o  Max  Criteria 


DFZE 

.  025°/hr 

DSZE 

. 025°/hr/g 

DIZE 

.03°/hr/g 

KIYE 

100  wg/g 

KIYHE 

ioo  wg/g 

KIXE 

loo  wg/g 

KIXHE 

ioo  wg/g 

DOZE 

. 02°/hr/g 

DFXE 

. 025°/hr 

KIZE 

ioo  wg/g 

KIZHE 

ioo  wg/g 

DFYE 

. 025°/hr 

DSXE 

. 025°/hr/g 

DSYE 

. 025°/hr/g 

DIXE 

. 03°/hr/g 

DIYE 

. 03°/hr/g 

DELTYX 

.00008  rad 

DELTZY 

.00020  rad 

KOZE 

300  wg 

KSZE 

50  wg/g 

KOZHE 

300  wg 

KSZHE 

so  wg/g 

KOXE 

ioo  wg 

KSXE 

so  wg/g 

KOXHE 

ioo  wg 

KSXHE 

so  wg/g 

DELTZX 

.00015  rad 

KOYE 

ioo  wg 

KSYE 

so  wg/g 

KOYIJE 

ioo  wg 

KSYHE 

so  wg/g 

KTXE 

. 0002°/hr  per 

KTYE 

•  0002Vhr  per 

KTZE 

•0002°/hr  per 

. 075 °/hr 
. 075°/hr/g 
. 09°/hr/g 
300  wg/g 
300  wg/g 
300  wg/g 
300  wg/g 
.06°/hr/g 
.075°/hr 
300  vg/g 
300  wg/g 
. 075°/hr 
. 075°/hr/g 
. 075°/hr/g 
. 09®/hr/g 
. 09*/hr/g 
.00025  rad 
.00060  rad 
900  wg 
150  wg/g 
900  wg 
150  wg/g 
300  wg 
150  wg/g 
300  wg 
150  wg/g 
.00045  rad 
300  wg 
150  wg/g 
300  wg 
150  wg/g 

/hr  . 0006°/hr  per  °/hr 

/hr  . 0006°/hr  per  ®/hr 

/hr  . 0006°/hr  per  °/hr 
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TABLE  3.  NUMBER  OF  IMU  CALIBRATIONS  AND  TIME  FRAMES 


IMU  S/N 

NO.  CAL. 

NO.  MONTHS 

TIME  FRAME 

REMARKS 

001 

26 

13.6 

3/76  -  5/77 

002 

9 

3.0 

2/76  -  5/76 

002 

18 

7.1 

5/76  -  12/76 

003 

6 

1.9 

9/76  -  11/76 

003 

9 

6.4 

1/77  -  7/77 

004 

7 

5.0 

9/76  -  2/77 

005 

10 

11.4 

6/77  -  6/78 

006 

3 

4.8 

5/77  -  10/78 

006 

4 

3.0 

12/77  -  3/78 

AD  FLT  4 

007 

5 

8.7 

6/77  -  3/78 

008 

3 

5.3 

6/77  -  11/77 

AD  FLT  2 

009 

1 

1.4 

9/77  -  11/77 

AD  FLT  3 

010 

1 

1.6 

8/77  -  9/77 

AD  FLT  1 

Oil 

2 

03 

. 

o 

7/77  -  11/77 

012 

1 

6.5 

8/77  -  2/78 

012 

5 

3.7 

4/78  -  8/78 

013 

7 

2.5 

2/78  -  5/78 

AD  FLT  5 

TOTALS 

117 

89.7 
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TABLE  4. 


DATE 

LOCATION 

ELAESED 

TIME 

(MONTHS) 

1 

DFZE  | 

DSZE 

DIZE 

KIYE 

06/28/77 

SKD 

0 (REF) 

08/05/77 

MMC 

1.27 

-.81 

1.93 

-.73 

-.35 

11/08/77 

RSA 

4.43 

.58 

1.56 

-1.53 

-.58 

11/29/77 

RSA 

5.13 

.  50 

1.48 

-.30 

-.45 

12/07/77 

RSA 

5.40 

-.11 

1.78 

-1.03 

-.43 

12/16/77 

RSA 

5.70 

0 

1.63 

-2.43 

-.83 

01/03/78 

RSA 

6.30 

.81 

1.38 

-.56 

-1.42 

01/11/78 

MMC 

6.56 

-.24 

1.62 

-1.45 

-1.10 

02/27/78 

MMC 

8.13 

.11 

1.28 

-1.51 

-1.55; 

03/23/78 

MMC 

1 

8.93 

-.02 

1.23 

-1.19 

-1.45 

06/06/78 

MMC 

11.43 

1.50 

2.70 

-9.70 

- i 

-1.8<j 

07T6cedxM$ 


9am#*  i 
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TABLE  4.  IMU  S/N  005  NORMALIZED  DATA 
REFERENCE  CALIBRATION: 

JUNE  28,  1977 


1 

DIZE 

KIYE 

KIYHE 

KIXE 

KIXHE 

DOZE 

DFXE 

KIZE 

KIZHE 

DFYE 

DSXE 

DSYE 

1 

■ 

3 

-.73 

-.35 

-.54 

-.06 

-.38 

-.70 

.08 

-.12 

.22 

-.58 

.16 

-.16 

3 

-1.53 

-.58 

-1.10 

n 

-1.28 

-1.00 

-.40 

-1.14 

.69 

-.74 

-.36 

.36 

i — 

-.30 

-.45 

-.87 

-.56 

-1.06 

-.50 

-.40 

1.08 

.64 

n 

.44 

5 

(8 

-1.03 

n 

-.78 

-.14 

-1.02 

-.20 

-.79 

1.12 

.64 

.39 

-.48 

_ 

.48 

-2.43 

-.83 

-1.15 

-.65 

-1.31 

.20 

-.56 

.76 

.39 

-.07 

-.48 

.48 

9 

-.56 

-1.42 

-1.84 

-1.51 

-2.02 

-1.60 

-.38 

— 

.65 

.29 

.63 

-.36 

.36 

m 

-1.10 

-2.24 

-.78 

1.55 

n 

-.36 

-.46 

-.72 

.10 

.12 

-.12 

8 

SO 

-1.55 

-1.86 

-1.28 

-1.72 

.  02 

-.72 

-.36 

-.55 

1.14 

-.28 

.28 

n 

-1.45 

-2.70 

-1.18 

1.12 

1.12 

-.74 

-.37 

-.62 

1.56 

-.36 

.36 

X 

SB 

-1.80 

-3.00 

-1.80 

.80 

.56 

-.09 

-.36 

-.56 

1.80 

.15 

-.15 

DATE 

LOCATION 

ELAPSED 
TIME 
(MONTHS ) 

DIXE 

. 

DIYE 

AYX 

AYZ 

IV 

■ 

06/28/77 

SKD 

0 (REF) 

■ 

08/05/77 

MMC 

1.27 

.10 

.10 

.07 

.01 

■B 

B§ 

11/08/77 

RSA 

mm 

-4.13 

-4.13 

.84 

.03 

mm 

m 

11/29/77 

RSA 

5.13 

-1.70 

-1.70 

1.16 

0 

■ 

12/07/77 

RSA 

5.40 

-4.60 

-4.60 

.68 

.01 

B 

12/16/77 

RSA 

5.70 

2.80 

2.80 

1.20 

0 

J 

B 

01/03/78 

RSA 

6.30 

-2.67 

-2.67 

.14 

.02 

B 

01/11/78 

MMC 

6.56 

-.07 

-.07 

.05 

0 

B 

02/27/78 

MMC 

8.13 

-.06 

— 

-.06 

.10 

.02 

B 

03/23/78 

MMC 

8.93 

-1.37 

-1.37 

.  34 

.03 

B 

06/06/78 

MMC' 

11.43 

-1.04 

-1.04 

.29 

.28 

B 
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TABLE  4.  (Continued) 


. 

AYZ 

KOZE 

KSZE 

KOZHE 

— 

KSZHE 

KOXE 

KSXE 

KOXHE 

KSXHE 

AZX 

KOYE 

.01 

-.31 

-1.47 

-.09 

.11 

-.34 

-.60 

-.06 

.14 

.41 

.56 

.03 

-.60 

-1.92 

-.25 

.69 

-1.45 

-3.64 

-.26 

-.72 

.54 

.65 

0 

-.07 

.80 

-.26 

.16 

-.08 

-.42 

J 

-.30 

-.63 

1.36 

L _ 

-.62 

.01 

-.19 

.29 

-.19 

.33 

.28 

.  09 

■  -1 

-.21 

-.59 

1.19 

-.47 

0 

-.06 

.22 

-.15 

.03 

-.10 

1 

-.11 

-.17 

-.57 

1.20 

-.57 

.02 

-.06 

.79 

-.23 

.15 

-.01 

-.32 

-.35 

-1.08 

- 

1.24 

-.35 

0 

-.32 

— 

-1.13 

-.21 

— 

-.23 

-.77 

-1.72 

.02 

-.02 

- 

.61 

.70 

.02 

-.19 

.21 

-.28 

-.04 

-.61 

-1.07 

.01 

-.33 

.68 

.64 

.03 

-.20 

.59 

-.36 

-.15 

-.88 

-1.75 

-.08 

_ _ _ _ 

-.19 

■  ■ 

.  66 

.56 

.28 

-.60 

-.60 

-.50 

.12 

-.60 

-1.30 

-.14 

-.34 

.70 

.40 

.  ■•***••  v- .  •  •  -  -••*-  ..  j 


TABLE  4.  (Concluded) 


DATE 

LOCATION 

ELAPSED 

TIME 

(MONTHS) 

KSYE 

KOYHE 

KSYHE 

06/28/77 

SKD 

0 (REF) 

08/05/77 

MMC 

1.27 

1.20 

-.10 

.26 

11/08/77 

RSA 

4.43 

1.27 

-.09 

.22 

11/29/77 

RSA 

5.13 

-  .  91 

.  16 

.71 

12/07/77 

- -  • 

RSA 

5.40 

— 

-.94 

■"  -  — 

.15 

.52 

12/16/77 

RSA 

5.70 

o 

r — 1 

1 

-.04 

.11 

01/03/78 

RSA 

6.30 

-.45 

-.18 

.22 

01/11/78 

MMC 

— 

6.56 

— 

1.71 

-.05 

.08 

02/27/78 

>*yc 

8.13 

1.00 

.10 

-.03 

03/23/78 

Vf  y.  Q 

8.93 

1.09 

.08 

.03 

06/06/78 

MMC 

11.43 

1.05 

.11 

.39 
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TABLE  4. 


(Concluded) 


LPSED 

I  ME 
KTHS) 

KSYE 

KOYHE 

KSYHE 

KTXE 

ktye  KTZE 

[REF) 

..27 

1.20 

-.10 

.  26 

-.70 

-.46  1.07 

1.43 

1.27 

-  .09 

.  22 

-.58 

.44  2.52 

- - 

5.13 

-.91 

.  16 

.71 

-.34 

-.74  1.95 

L  40 

-.94 

.15 

.  52 

.53 

.10  2.68 

i.70 

-1.09 

-.04 

.11 

.92 

-.65  2.30 

[.30 

-  .45 

-.18 

.22 

1.04 

4.19  3.10 

i.56 

1.71 

-.05 

.08 

.16 

-.34  2.86 

1.13 

1.00 

.10 

-  .03 

1.80 

-.10  3.10 

.93 

1.09 

.08 

.03 

1.94 

.06  3.42 

.43 

1.05 

.11 

.39 

1.80 

.31  3.10 
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RATION  DATA  -  PARAMETER  RANGE  (INDIVIDUAL  IMUs ) 
RE  FREQUENCY  OF  OCCURRENCE 


KIYE 

SIGMA  RANGE 


KIYHE 

SIGMA  RANGE 


KIXE 

SIGMA  RANGE 


KIXHE 

SIGMA  RANGE 


TABLE  5.  (Continued) 


,«**#.  Ur 


4*3, . .. 


— 

IMU 

S/N 

DOZE 

SIGMA  RANGE 

DFXE 

SIGMA  RANGE 

KIZE 

SIGMA  RANGE 

KIZHE 

SIGMA  RANG 

E 

+] 

+  1  |  +2 

+  3 

>  +  3 

+  1 

a 

a 

ggl 

a 

a 

a 

>  +  3 

a 

m 

+  3 

>  +  3 

001 

— 

22 

3 

0 

1 

10 

2 

9 

a 

a 

a 

0 

22  ■ 

2 

2  ! 

0  1 

19 

002 

25 

2 

o 

0 

27 

0 

0 

0 

11 

n 

2 

3 

27  i 

o  ! 

°  i 

i 

0  1 
1 

27 

003 

15 

0 

0 

0 

7 

1 

4 

3 

14 

i 

0 

0 

14 1 

1 

0  i 

i  i 

° 

8 

004 

7 

0 

0 

0 

7 

0 

o 

0 

7 

a 

0 

0 

7 

0 

0  i 

1 

o  1 

4 

005 

8 

2 

0 

0 

10 

0 

0 

0 

7 

3 

0 

0 

m 

0 

0  : 

0 

7 

006 

3 

3 

1 

0 

7 

0 

o 

o 

4 

2 

1 

0 

7 

0 

0  ! 

0 

7 

007 

5 

0 

0 

0 

0 

2 

1 

2 

5 

0 

0 

0 

1 

2 

0 

2 

1 

008 

2 

1 

0 

0 

1 

0 

2 

0 

3 

0 

0 

0 

3 

0 

0 

0 

3 

009 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0  : 

0 

0 

1 

010 

0 

1 

0 

0 

1 

0 

0 

o 

1 

0 

0 

0 

1 

o 

i 

0 

0 

1 

Oil 

2 

0 

0 

0 

0 

0 

1 

1 

2 

0 

0 

0 

2 

0 

0 

0 

0 

012 

6 

0 

0 

0 

6 

0 

0 

0 

5 

1 

0 

0 

4 

2 

0 

0 

6 

013 

6 

1 

0 

0 

2 

3 

2 

0 

2 

2 

3 

0 

7 

0 

_ 

0 

0 

— 

7 

TOTALS 

101 

14 

1 

1 

79 

11 

12 

15 

6 

3 

3 

2 

91 

% 

86 

12 

1 

1 

68 

9 

10 

13 

74 

18 

5 

3 

90 

5 

3 

2 

78 
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TABLE  5.  (Continued) 
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TABLE  5.  (Continued) 


DELTYX 


DELTZY 


KOZE 


KSZE 


KOZHE 


RA 

NGE 

3IGMA  RANGE 

3IGMJ 

\  RANGE 

STOM 

A  RANGE 

SIGMA  RAI 

+  3 

>  +3 

+1 

+  2 

+  3 

>  +3 

i1 

t2 

+  3 

>  +3 

+1 

+  2 

+  3 

>  +3 

+  1 

+  2 

+  3 

0 

0 

26 

0 

0 

0 

26 

0 

0 

0 

7 

7 

12 

0 

26 

0 

0 

0 

0 

27 

0 

0 

0 

27 

0 

0 

0 

24 

3 

0 

0 

27 

0 

0 

0 

0 

15 

0 

0 

0 

14 

0 

0 

12 

3 

0 

0 

15 

o 

0 

0 

0 

7 

0 

0 

0 

7 

0 

0 

0 

1 

5 

1 

0 

7 

0 

0 

0 

0 

10 

0 

0 

0 

10 

0 

0 

0 

7 

3 

0 

0 

10 

0 

0 

0 

0 

7 

0 

0 

0 

7 

0 

0 

0 

3 

4 

0 

0 

7 

0 

0 

0 

0 

5 

0 

0 

0 

5 

0 

0 

0 

4 

1 

0 

0 

3 

1 

0 

0 

0 

3 

0 

0 

0 

1 

2 

0 

0 

i 

0 

0 

2 
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One  method  that  is  widely  used  to  examine  data  for 
trending  effects  is  linear  regression. 2  a  straight  line 
of  the  form 

Y  =  A  +  BX  (1) 


is  fitted  to  the  data  such  that  the  sum  of  squares  of  ver¬ 
tical  deviations  of  the  observed  data  from  this  line  is 
smaller  than  the  corresponding  sum  for  any  other  straight 
line.  The  regression  line  is  therefore  a  least-square 
estimate  of  the  unknown  true  line.  Equation  (1)  is  the 
familiar  slope/intercept  form  with 


NEXY-EXIY 
NEX2-  (EX)  2 


and 


(2) 


! 


ZY-BZX 

N 


(3) 


in  which 

N  *  number  of  data  point  pairs 
X  ■  independent  variable 

Y  =  dependent  variable  (observed  data,  each  value 
paired  with  a  corresponding  value  of  X) 

2  =  summation  of  indicated  variables  over  all  data 
points  from  1  to  N. 

After  determination  of  A  and  B  it  is  a  simple  matter 
to  use  Equation  (1)  to  provide  an  estimate  of  Y  for  a  n6w 
specified  value  of  X.  Alternately,  X  can  be  calculated  for 
any  value  of  Y: 

X  =  ^  (4) 

T.  US  Nival  Ordnance  Test  Station;  Statistics  Manual  by 
Crowe,  E.  L. ,  David,  F.  A.,  and  Maxfield,  Margaret 
W.,  China  Lake,  California,  1955,  Chapter  6. 


frecediA'G  ^ 
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Equation  (4)  provides  an  optimistic  estimate  of  X 
because  the  actual  Y  data  is  scattered  about  the  regression 
line.  The  scatter  is  measured  by 


e2  - 
“yx  "  N-2 


(sJ-b2s*  ) 


(5) 


and  the  positive  square  root  is  sometimes  called  the 
standard  error  of  the  estimate. 

In  Equation  (5) 


q2  _  NEY2- (£Y) 2 
y  N(N-l) 


and 


_2  _  NZX2-  (IX)  2 
x  N(N-l) 


<«) 


(7) 


which  are  the  Y  and  X  variances  about  their  respective 
means. 


If  the  Y  data  corresponding  to  the  various  Xs  is 
normally  distributed  about  the  true  regression  line,  two 
lines  drawn  parallel  to  the  calculated  regression  line  at 
vertical  distances  +‘  SyX  envelop  68  percent  of  the  data. 
These  lines  are  defTnea  by 


Y  =  (A+S  ) +BX  (8) 

and 


Y  =  (A“Syx)+BX  (9) 

and  are  illustrated  in  Figure  A-l.  If  it  is  desired  to 
estimate  a  maximum  X  for  an  assumed  value  of  Y,  e.g.,  Y^ 
in  Figure  A-l,  the  conservative  approach  would  be  to  use 
Equation  (8)  or  the  top  line  of  Figure  A-l  to  account  for 
some  of  the  variation  about  the  calculated  regression  line. 
Thus  X^  would  be  the  conservative  estimate. 


Ml 
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Further  conservatism  may  be  applied  to  the  estimate 
by  consideration  of  the  standard  error  of  the  slope  which 
is  given  by 


c  =  _ _ 

B  SX/N^I 


Equation  (10)  allows  for  variations  of  the  slope  B 
and  is  used  to  define  two  lines  given  by 


Y  =  A+ (B+Sg) X 


Y  *  A+(B-Sb)X  (IS 

which  are  illustrated  in  Figure  A-2.  Again  Xi  would  be 
the  conservative  estimate- ?or—a"'gTven  value  of  Yj_. 

Allowance  for  variation  of  slope  and  variation  of 
data  about  the  line  of  chosen  slope  can  be  made  by  combin¬ 
ing  elements  of  Equations  (8)  and  (11)  and  Equations  (9) 
and  (12)  to  give 


Y  =  (A+Syx)  +  (B+Sb)X 


Y  =  (A-Syx)  +  (B-Sb)X 


Plots  of  Equations  (13)  and  (14) ,  as  shown  in  Figure 
A-3 ,  provide  a  funnel  about  the  original  regression  line 
and  should  enclose  practically  all  of  the  observed  data. 

Equation  (13)  provides  the  most  conservative  estimate 
of  X  for  any  of  the  cases  considered  when  the  slope  B  is 
positive  as  indicated  in  Figure  A-3. 

Similarly,  when  the  slope  is  negative.  Equation  (14) 
provides  the  conservative  estimate. 


The  criterion  for  the  estimation  of  a  conservative 
value  of  X  when  Y  *  +Yj  is  thus  established  as 


y  ±Yr  <*±Syx> 
1  =  (B+Sg) 


(15) 


in  which  the  sign  of  Y^,  SyX  and  S3  is  taken  the  same 
as  the  sign  of  the  slope  B. 

In  the  case  of  the  PII  IMU  calibration  trend 
analysis,  an  estimate  of  X  (months)  is  sought  which  indi¬ 
cates  when  +Y2  =  +3  sigma  since  the  +3  sigma  limits 
specify  the  unqualified  requirement  7or  a  new  calibration 
under  AD  ground  rules. 

The  linear  regression  line  as  given  in  Equation  (1) 
was  established  for  each  parameter  from  all  IMUs  having 
six  or  more  calibrations  over  a  period  of  six  or  more 
months  (serial  numbers  001,  002,  003  and  005).  Equation 
(15)  was  used  to  calculate  a  conservative  estimate  of  X  * 

Xi  (months)  which  would  be  expected  to  elapse  before  a 
parameter  exceeded  its  +3  sigma  limits  thus  giving  an 
indication  of  the  amount  of  time  that  may  be  allowed 
between  calibrations. 

Figures  A- 4  through  A- 3 7  show  data  obtained  for  each 
parameter  from  IMU  serial  number  005  and  the  corresponding 
regression  lines  given  by  Equation  (1).  Also  shown  are  the 
lines  specified  by  Equations  (13)  and  (14). 

Other  statistics  computed  from  calibration  data  from 
those  IMUs  analyzed  are: 

(1)  Correlation  coefficient: 

R  -  BSx/Sy  (16) 

with  B  given  by  Equation  (2)  and  Sy  and  Sx  given  by  the 
respective  square  roots  of  Equations  (6)  and  (7). 

(2)  Root  mean  square  error  of  the  dependent  variable: 

Y  ilS)  =  ( IY2/N)  ^  (17) 
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with  symbol  definitions  the  same  as  those  given  for 
Equations  (2)  and  (3). 

(3)  Average  error  of  the  dependent  variable: 


Y(AVG)  *  EY/N  (18 


with  symbol  definitions  the  same  as  those  given  for 
Equations  (2)  and  (3). 

Computed  values  of  all  statistics  for  each  calibra¬ 
tion  parameter  from  IMUs  with  serial  numbers  001,  002,  003 
and  005  along  with  conservative  estimate,  X,  and  estimate, 
Xz,  from  the  original  regression  line,  are  compiled  in 
Tables  A— 1  through  A-4  respectively. 

Resulting  values  of  the  conservative  estimate,  Xj., 
are  summarized  in  Table  A-5  which  also  shows  the  values  of 
Xi  for  each  parameter  averaged  over  the  four  IMUs 
analyzed. 


Linear  regression  lines  showing  +  S  bounds. 


Figure  A-2.  Linear  regression  lines  showing  +SB  bounds. 
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Figure  A-10.  KIXHE 
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